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In the supersymmetric framework, a higgsino asymmetry exists in the universe before the electroweak 
phase transition. We investigate whether the higgsino is a viable asymmetric dark matter candidate. 
We find that this is indeed possible. The gauginos, squarks and sleptons must all be very heavy, such 
that the only electroweak-scale superpartners are the higgsinos. The temperature of the electroweak 
phase transition must be in the (1 — 10) GeV range. 



Introduction. The matter content of our universe is 
made of two main components: the dark matter (DM), 
with Oqm ~ 0.20, and baryons, with fib ~ 0.04. Intrigu- 
ingly, neither of these numbers can be explained within 
the Standard Model (SM) of particle physics. The most 
intensively studied scenarios involve very different mech- 
anisms to explain these two numbers. The DM relic 
abundance is explained by a freeze-out of a weakly in- 
teracting massive particle number density that occurs 
when its annihilation rate becomes slower than the ex- 
pansion rate of the universe. The baryon relic abundance 
is explained by an asymmetry between baryons and anti- 
baryons. Under these circumstances, there is no explana- 
tion of the fact that the energy densities of the DM and 
the baryons are surprisingly close to each other, which is 
then just a coincidence. It would be more satisfying if it 
could be naturally explained. This can be the case if the 
DM density were also the result of an asymmetry, rather 
than of freeze-out [lTfT2]. This type of scenarios comes 
under the name of asymmetric dark matter. 

One of the best motivated extensions of the SM is the 
Minimal Supersymmetric Standard Model (MSSM). Ex- 
tending this framework, there are several ways to gener- 
ate the baryon asymmetry. Perhaps the most plausible 
one is via leptogenesis (for a review, see [13] )• Prior to the 
Electro Weak Phase Transition (EWPT), when the hig- 
gsino is of Dirac nature, leptogenesis generates also Higgs 
and higgsino asymmetries that are initially of similar size 
to the lepton asymmetry (see e.g. Q3]). Regardless of 
the source of the baryon and higgsino asymmetries, how- 
ever, the conditions of chemical equilibrium imply that, 
at temperatures above the EWPT, a non-zero baryon 
asymmetry requires that there is also a non-zero higgsino 
asymmetry. Could such a higgsino asymmetry survive 
and lead to asymmetric higgsino dark matter after the 
EWPT, when the higgsino becomes of Majorana nature? 
This is the question that we address in this work. 

Imagine that the following set of conditions applies: 
1) the higgsino is the lightest supersymmetric particle; 



2) before the EWPT, higgsino number-changing inter- 
actions become slow enough that the higgsino asym- 
metry remains constant; 3) after the EWPT, higgsino- 
antihiggsino oscillations are slow or, if they are fast, the 
higgsino-antihiggsino annihilation rate is slow. Then, a 
rather large relic higgsino density can survive. In this 
work, we study the constraints on the supersymmetric 
spectrum and on cosmological parameters that follow 
from imposing this set of conditions. 
The higgsino asymmetry. The asymmetry in a par- 
ticle x number density, An x = (n x — n x ), is related to its 
chemical potential, fi x , via (for fi x <C T) 



An7 = 9*r^ K[zxl 



(i) 



where g x is the number of internal degrees of freedom of 
the particle /T, K(z x <1) = 2 (1) for bosons 

(fermions) and K(z x 3> 1) is exponentially suppressed. 
We are interested in relating the higgsino asymmetry to 
the baryon asymmetry. Under the conditions that will 
be of interest to us (see below), all of the sfermions are 
much heavier than the higgsinos and the sfermion number 
densities are negligible because the relevant K (2) factors 
are exponentially suppressed. On the other hand, for the 
quarks and leptons K(z) — 1. With these assumptions, 
the baryon and lepton asymmetries are given by 

AFb = — (2^ Q +M„ + Md), AF L = — (2Mt+A*e), (2) 

where ^ = J2i A^Vi (i = lj 2, 3 is a generation index). We 
define the comoving asymmetry via AY X = An x /s, where 
s is the entropy density. We follow the derivation of Rcfs. 
[HI [16]. Imposing the conditions of fast gauge, Yukawa, 
and sphaleron interactions, the Majorana nature of the 
gauginos, the Dirac nature of the higgsinos (prior to the 
EWPT), and hypercharge neutrality, we obtain 
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2K(h) 



12 + 3[K(h u ) + K{h d ) + 2K(h)} 



(3) 
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Eq. ([3]) relies on super-equilibrium (SE), that is chem- 
ical equilibrium between Higgs particles and higgsinos. 
At the end of the epoch of SE, a higgsino asymmetry is 
conserved until the EWPT. 

As long as the higgsinos arc relativistic, the right hand 
side of Eq. ^ is (in absolute value) in the range (0.10 — 
0.15). If the higgsinos become non-relativistic while SE 
persists, then AY^ is quenched by the K(h) factor in 
Eq. ([3]). Such quenching of AY^ is not allowed if higgsinos 
are to provide the DM. Thus, we are led to impose that 
SE must be broken while higgsinos are still relativistic, 
freezing AYr at a value 

ay^-lq- 11 . (4) 

In the next section we find the conditions on the super- 
symmetric spectrum that would lead to early breakdown 
of SE. 

Non-Super-Equilibrium (NSE). Our goal is to find 
the conditions on the particle spectrum such that NSE 
occurs while the higgsinos are relativistic, Tnse > M- We 
use r < H as the rough criterion for an interaction to 
be out of equilibrium, where the Hubble expansion rate 
is given by H « 10T 2 /mpi and the interaction rate is 
given by r = n(av), n being the number density of target 
particles. 

Higgsino number is violated in Higgs-Higgs scattering 
(hh — > hh) and Higgs-antihiggsino scattering (hh c — > 
h c h). These processes arise from the effective Lagrangian 

- £ c ff = -^-huh u h* u h* u + ^-h d h d h* d h* d , (5) 
A« A d 

generated at tree level by gaugino and at one loop by 
quark-squark diagrams. Requiring that h u h u — > h u h u is 
not in equilibrium at T ~ /i puts a lower bound on A„, 



k u > 3 x 10 9 GeV ( -^ J —) 
VlTeV/ 



1/2 



(6) 



If h d is not heavy and decoupled, then a similar bound 
holds for Ad,. 

The gaugino contributions to ^ are given by 

1 _ 1 _ g' 2 g 2 
A^~T d ~Mf 1 + 8M2' 

The stop and sbottom contributions are given by 



(7) 
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Eq. (JsJ) , therefore, implies 
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M, > 10 s GeV 
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(9) 



Finally, decoupling the decays and inverse decays of 
sfermions to higgsino-fermion [16] at Tnse > M requires 
all of the sfermions to be heavy, m;> (10 — 40) /i. Here, 
the stronger bound refers to top squarks and the weaker 
bound refers to the electron super-partners. 

We note that since, as shown in the next section, 
A* > Tewpt, the above spectrum guarantees that SE ends 
before the EWPT. 

The spectrum. The higgsino spectrum includes two 
neutral mass eigenstates, with mass splitting Atoo |18j . 
and a charged higgsino, split by Am + from the lightest 
neutral higgsino [TH1 [20] . As concerns Am+ , this quan- 
tity does not violate higgsino number and in our sce- 
nario it is therefore dominated by gauge loops, Am + w 
^aws^ymz ~ 350 MeV. The results of the previous 
section have, however, interesting implications for Amo- 
The operators of Eq. (|5| lead to 

Am = A„ + Ad, (10) 
A» = (h u ) 2 /A u , A d = (h d ) 2 /A d . 

Eq. (|6]) implies A„ < 10 keV. An analogous bound, 
A d < (10 keV)/ tan 2 j3, applies unless the second Higgs 
doublet, h dl is much heavier than the higgsinos. We now 
show that experimental constraints lead us to choose the 
second possibility, namely mh d S> [i. 

Constraints on the higgsino spectrum in our scenario 
come from direct and indirect DM searches. The dom- 
inant higgsino-nucleus interactions are spin-independent 
inelastic interactions (see, e.g., [20-22 ). The cross sec- 
tion with a single nucleon, in the limit of zero mass split- 
ting, is a n w 10" 38 (10~ 36 ) cm 2 for fj, = 100 (1000) GeV. 
This is orders of magnitude above current bounds. The 
only way to evade these bounds (and maintain higgsino 
LSP) is by a large enough mass splitting that will make 
the inelastic scattering kinematically forbidden [2"3Tl25| . 
The minimum mass splitting required is a function of the 
higgsino mass, but in the entire range of interest for fi 
it is smaller than 400 keV. We learn that the sbottom 
contribution to Amo must be substantial, A^ > 400 keV. 
As a result, the only configuration consistent with early 
breakdown of SE, Tnse > Mi is one with a very heavy 
second doublet, my id ^> \i. 

As can be seen from Eqs. 



Am > 1 MeV 
constrains the sbottom sector to satisfy 

sin20 fc >2xlO- 2 (^\ 



and (10), requiring 



(11) 



(12) 



Finally, with Amo S> 10 eV, the present higgsino popu- 
lation is symmetric (see below), and therefore annihilates 
and may provide signals in indirect searches for DM. The 
annihilation cross section into WW and ZZ pairs is given 
by [2£] 
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(13) 
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The bound on the cross section from the Fermi-LAT data 
[27j . when compared with Eq. (13), requires 



(i > 190 GeV. 



(14) 



Oscillations, damping and expansion. At the 

EWPT, the Higgs acquires a VEV, and the higgsinos mix 
with the gauginos. The resulting propagation eigenstates 
change from Dirac to Majorana fermions, and oscillations 
begin. On the other hand, incoherent interactions with 
the plasma continue and damp the oscillations. In this 
section, we study the time evolution of the higgsino sys- 
tem (h, h) under the simultaneous effects of oscillations, 
annihilations, damping and the expansion of the universe. 
To do so, we employ the formalism of the density matrix. 
This formalism was originally developed to study neutri- 
nos |28H31| and we adapt it to our case. Our equations 
are consistent with those of Ref. [35] which deals with 
closely related issues (see also |3H|). The quantum rate 
equations for Y" M = n^/s (where no ± n 3 are the number 
densities of the higgsinos and anti- higgsinos, and n\ =F irii 
are the off diagonal elements of the density matrix) read 



dlogz 




dlogz 
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where z = fi/T. The damping (or decoherence) factor 
D, and the effective matter potential V, are given by 



D 



v),V = 8((3)a w r]BT 3 /(irm 2 w ). 



f 



While D is proportional to the elastic scattering cross 
section and to the total number density of the massless 
fermions in the plasma, V is proportional to the elas- 
tic scattering amplitude and to the fermion-antifermion 
asymmetry. In the limit of large damping, the effective 
rate of oscillation is r osc ~ (Atuq) 2 /D, as obtained in 
[T2"] . The G factor measures the ratio between the an- 
nihilation cross section with and without including co- 
annihilations. 

The initial conditions, at the EWPT, are the following: 



Y = n- h /s, Y 1 =Y 2 = 0, *3 = -10- 



(15) 



where is the solution of the Boltzmann equations (at 
the EWPT) with constant asymmetry, and the value of 
Y 3 = AY^ is taken from Eq. Q. 

Asymmetry-assisted higgsino DM. We aim to solve 
for Yo(co) which gives the final total number density in 
higgsinos, and for 13(00) which gives the final higgsino 
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FIG. 1: The required Amo as a function of /U, with Tewpt 
100 GeV. 



asymmetry. In particular, to provide f^DM h 2 
scenario needs to fit 



0.11, 



-v^obs 
*0 



7.6 x 10 



-13 



1 TeV///) . 



(16) 



We consider three free parameters: //, Amo and Tewpt 
and ask whether there is a range of these parame- 
ters where such a fit is achieved. If the asymmetry is 
washed out before the symmetric decoupling tempera- 
ture, « /i/25, then the present DM relic density is 
the standard symmetric one, as if an initial asymmetry 
was never generated: 



F sym w 5.9 x 10 



-13 



(/VI TeV) . 



(17) 



Comparing Eqs. (16) and (17) we learn that for // ~ 1.1 



TeV, higgsinos can account for DM without asymmetry. 
Since in the presence of an asymmetry the total number 
density is always larger than in the symmetric case, this 
puts an upper bound of 1.1 TeV on the higgsino mass. 
Thus, the range of interest is 190 GeV < // < 1.1 TeV. 
Within this range, the lighter the higgsino, the larger the 



asymmetry that is required in order to satisfy ( 16 ) 
We can obtain Y (oo 



> F n & 



if either of the follow- 



ing two conditions applies: 1) Amo is small enough that 
the oscillations are slow, and at least part of the higgsino 
asymmetry survives down to T < T^ n ; 2) The EWPT 
occurs late enough that annihilations are already slow 
when oscillations begin, Tewpt < ^dc™- To study the 
first possibility, we fix T EW pt = 100 GeV, and solve for 
Fo(oo) as a function of /i and Amo- We find that the 
lighter the higgsinos are, the smaller the mass splitting 
that is required to provide the DM abundance. The rea- 
son is that smaller /i leads to smaller T^™, and the asym- 
metry is required to survive to later times. The required 
Am as a function of /i is shown in Fig. [T] It ranges be- 
tween (10~ 2 - 10~ 8 ) kcV for /x in the range (1000 - 200) 
GeV. Such a small mass splitting is excluded by direct 
DM searches. We conclude that it is impossible for the 
asymmetry to survive once the EWPT takes place. 

Second, we fix Amo = 1 MeV (the results are not 
sensitive to changes in Amo m the range where it is 
not excluded by direct searches), and solve for Fq(oo) 
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FIG. 2: The required Tewpt as a function of /i, with Amo = 
1 MeV. The approximate analytic solution ( 18 1 is shown in 
dashed line. 



as a function of \x and Tewpt- Our numerical result for 
the required Tewpt as a function of [i is shown as the 
smooth line in Fig. [2] Above the line, the final DM abun- 
dance is too low. Below the line, the final abundance 
is too high. The required temperature ranges between 
(30 - 0.1) GcV for fi in the range (1000 - 200) GeV. For 
T > 3 GeV(/i/l TeV) 2 , we can obtain an approximate 
analytic solution for the required Tewpt, 



Tewpt = 33 GeV 



1 TeV 



(18) 



shown as the dashed line in Fig. [2] 

We conclude that a viable scenario of asymmetric hig- 
gsino dark matter could have occurred as follows: a hig- 
gsino number asymmetry of a size that is about a factor of 
10 smaller than the baryon asymmetry exists before the 
EWPT. At the EWPT, the asymmetry is quickly washed 
out due to higgsino-antihiggsino oscillations. The result- 
ing symmetric higgsino population is (for masses below 
TeV) much larger than the would-be population without 
an initial asymmetry. It survives if the phase transition 
occurs at a temperature that is somewhat low, of or- 
der (1-10) GeV. We note that such a low temperature 
requires a rather strong phase transition, potentially ne- 
cessitating new degrees of freedom other than those of 
the MSSM. Further analysis of this requirement is be- 
yond the scope of the current paper and we postpone it 
to future work. 

Conclusions. Within the framework of the MSSM, we 
ask whether the higgsino could be a viable asymmetric 
dark matter candidate. We find that the answer is in 
the affirmative, provided that the following constraints 
on the supersymmetric spectrum are satisfied: 

• Electroweak gauginos are heavier than 10 8 GeV; 

• Sfermions are heavier than 10 4 GeV; 

• The stop mixing angle is small, and the sbottom 
mixing angle is large; 

• Higgsinos are in the range (200-1000) GeV. 



In addition, the temperature of the electroweak phase 
transition must be somewhat low, of order (1-10) GeV. 

The supersymmetric spectrum is somewhat reminis- 
cent of split supersymmetry models [3H |3S] . The super- 
symmetric flavor problem is solved. Grand Unification 
remains a viable possibility. Supersymmetry does not 
solve the fine tuning problem. It does however explain 
both the baryon asymmetry and the dark matter abun- 
dance, and relates the two. The initial source of both 
asymmetries could be leptogenesis. 

This scenario, where the only new particles at the elec- 
troweak scale are the higgsinos, poses a challenge to the 
LHC. Work on experimental and observational signals is 
in progress. 
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